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We briefly review the solution of three ensembles of non-Hermitian ran- 
dom matrices generalizing the Wishart-Laguerre (also called chiral) ensem- 
bles. These generalizations are realized as Gaussian two-matrix models, 
where the complex eigenvalues of the product of the two independent rect- 
angular matrices are sought, with the matrix elements of both matrices 
being either real, complex or quaternion real. We also present the more 
general case depending on a non-Hermiticity parameter, that allows us to 
interpolate between the corresponding three Hermitian Wishart ensembles 
with real eigenvalues and the maximally non-Hermitian case. All three 
symmetry classes are explicitly solved for finite matrix size iV x M for all 
complex eigenvalue correlations functions (and real or mixed correlations 
for real matrix elements). These are given in terms of the corresponding 
kernels built from orthogonal or skew-orthogonal Laguerre polynomials in 
the complex plane. We then present the corresponding three Bessel kernels 
in the complex plane in the microscopic large- iV scaling limit at the origin, 
both at weak and strong non-Hermiticity with M — iV > fixed. 

1. Introduction and motivation 

The Wishart ensemble was the first ensemble of random matrices, intro- 
duced in the study of randomized rectangular times series matrices C in or- 
der to study the spectral properties of its symmetric and positive covariance 
matrix C^C. Later the Wigner-Dyson ensembles were formulated as mod- 
els for randomized Hamiltonians H of heavy nuclei to explain some of the 
spectral properties of the matrix H which is Hermitian (or real symmetric 
or quaternion self dual). Its non-Hermitian generalizations were introduced 
and studied immediately after by Ginibre [2T], but not so for the Wishart 
ensembles which had to wait a few decades. 
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In order to construct such a generalization it is useful to compare to 
the Ginibre ensembles and to think of them as a two-matrix problem. By 
dropping the Hermiticity constraint an independent anti-Hermitian matrix 
A is added to the Hermitian matrix H, and the now complex eigenvalues of 
the sum H + A oi the two are studied. How can one repeat this for Wishart 
matrices, where the matrix C is already without symmetry? Again one can 
introduce a second independent random matrix D, and study this time the 
complex spectrum of the product DC, where D is no longer the transpose 
(or Hermitian conjugate) of the first matrix, ^ D, but has the same 
rectangular dimension. 

Interestingly the solution of this two-matrix problem did not grow out 
of statistical applications, where the problem has appeared in terms of so- 
called time lagged, asymmetric covariance matrices e.g. in \28\ [29l lllj . 
Non-Hermitian Wishart ensembles also appeared naturally as the described 
two-matrix problem in the study of the Dirac operator spectrum of Quan- 
tum Chromodynamics (QCD) with chemical potential /x. Here the ensemble 
with complex (/3 = 2) [32], quaternion real (/3 = 4) [2], and real (/3 = 1) 
[H E] matrix elements were first introduced and solved for finite and infinite 
matrix dimensions by Osborn, the author and his coworkers, respectively. 
In particular this includes the generalization of the Bessel kernels in the mi- 
croscopic origin scaling limit into the complex plane for all three ensembles. 
We refer to [38] for the most recent review on the topic of random matrix 
applications to QCD. Because the solution of these three ensembles can be 
expressed in terms of Laguerre polynomials in the complex plane, and be- 
cause these ensembles display chiral symmetry (see e.g. [38]) they are also 
called non-Hermitian (or complex) Laguerre or chiral ensembles. The link 
to statistical applications of non-Hermitian Wishart ensembles was reem- 
phasized more recently in [26]. Here and independently in [T2] the spectral 
density generalizing the Marchenko-Pastur distribution into the complex 
plane was computed. In [12] the product of having also more than two 
rectangular matrices was considered. 

The two-matrix models constructed and solved for applications to QCD 
are much more general, for two reasons. First, they depend on a non- 
Hermiticity parameter (the chemical potential //) that allows to smoothly 
deform the Hermitian Wishart (-Laguerre or chiral) ensembles into the max- 
imally non-Hermitian Wishart ensembles of two independent matrices de- 
scribed above. In the Ginibre ensembles these deformations also exist and 
are called elliptic, or Ginibre-Girko ensembles. The parameter fi allows to 
study deformations of the kernels known to be universal from the Hermi- 
tian setting. Second, more source terms were added to these parameter 
dependent two-matrix models, by inserting an arbitrary but fixed number 
of characteristic polynomials (as additional determinants of the Dirac oper- 
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ator) into the measure. All complex eigenvalue correlation functions were 
computed in this more general setting in [32; 1] for /3 = 2, in [2l|3] for /3 = 4 
and in |8l |9] for /3 = 1. We also mention that the /3 = 1 symmetry class 
at maximal non-Hermiticity appears in the superconducting phase of QCD 
with two colors [24J. In this short presentation we will focus on the first 
aspect, the dependence on the non-Hermiticity parameter and recall the full 
solution in terms of (skew) orthogonal polynomials in the complex plane. 
The additional characteristic polynomials can then be easily implemented 
by modifying the corresponding kernels (see e.g. in [HIS])- 

Finally we would like to mention some related developments in non- 
Hermitian random matrices. First of all non-Hermitian generalizations of 
Wishart ensembles were first considered as one-matrix models [37t [23] where 
the non-Hermiticity is provided by a constant matrix shift. These models 
are much more difficult to be solved in general, and only the microscopic 
density for f3 = 2 was determined in |36j . 

Instead of considering the complex eigenvalues of the product of two 
matrices, in [16] those of the ratio of two matrices were studied, which re- 
late to a Cauchy distribution. When generalizing from one to two Wishart 
matrices, one can also consider the positive Hermitian (real symmetric) 
combination (DC)^ DC , as was done in for instance in [33], including more 
matrices. Again this generalizes the Marchenko-Pastur density, this times 
for real eigenvalues, finding generating functions relevant in combinatorics. 
More generally speaking, spectral properties of the product of quadratic ran- 
dom matrices - Hermitian or non-Hermitian - have been studied by several 
authors in the literature, and we refer to [13] as well as to the contribution 
[13] to these proceedings and references therein. The three Ginibre ensem- 
bles, and the three non-Hermitian Wishart ensembles reviewed here are not 
the only possible non-Hermitian random matrices one can consider. For an 
ordering principle we refer to [TUl EH] which include these 6 classes out of 33 
non-Hermitian ones. For recent reviews on non-Hermitian random matrices 
we refer to [271 [32] ■ 

This review is organized as follows. In Sect. [2] we recall the definition 
of generalized non-Hermitian Wishart ensembles, discuss their relation to 
standard Wishart ensembles and give their complex eigenvalues representa- 
tions. Sect. [3] gives a list of the three sets of (skew) orthogonal Laguerre 
polynomials in the complex plane and their kernels that allows to compute 
all complex (and real) eigenvalue correlation functions. The large- limit is 
sketched in Sect. [H starting with the elliptic law for the Dirac eigenvalues 
and then listing the three Bessel kernels in the microscopic origin limit at 
weak and strong non-Hermiticity. A short discussion on universality follows 
in Sect. [3 



4 



Allv2 PRINTED ON JANUARY 20, 2013 



2. Complex eigenvalue representation of the partition function 

We begin by recalling the definition of the standard Wishart ensembles 
we wish to generalize. Its partition function is given by 



dC exp 



--Tr(CtC) 



(2.1) 



where C is a rectangular matrix of size N x (A^ + v) with real (/3 = 1), 
complex (/3 = 2) or quaternion real (/? = 4) matrix element without further 
symmetry. When going to a real eigenvalue basis for these ensembles there 
are two choice^]. Either we can consider the N positive eigenvalues Aj of 
the positive Hermitian matrix C^C, which we will call Wishart eigenvalues 
(or Wishart picture). Or we can consider the 2N + v eigenvalues Xi of 

the Dirac matrix p = { ^Jl „ ^ ] which we call Dirac eigenvalues (or 

Dirac picture). When comparing the two characteristic equations, 



N N 

det[A - C^C] = JJ(A - \i) , vs. det[A -p\ = \'' JJC^^ " ) ' (2-2) 

i=l i=l 

it is clear that we have v zero- and 2N non-zero Dirac eigenvalues that come 
in ± pairs. The simple substitution Aj = xj will lead us from one picture 
to the other for the non-zero eigenvalues. When computing the standard 
Jacobian for the diagonalization we obtain the partition function in terms 
of eigenvalues as follows 



OO ^ R 



\[dX.Xr^'^''e~^'^ |A^({A,})|^ 

i=l 
poo N 

/ lldx,\x.f('^+'^-'e-^^^ |A^({x|})|'^. (2.3) 



It is given in both the Wishart and Dirac picture, and we have defined the 
Vandermonde determinant as 



Aiv({A,})= n (Afc-AO 

N>k>l>l 



det 

l<kA<N 



(2.4) 



^ We note in passing that on the level of matrix elements eq. H2.1|l is identical to the 
Ginibre ensembles for u = Q. The difference is that there the complex eigenvalues of 
the matrix C are studied. 
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Wishart picture 


Dirac picture 


a) 


Vj > real 


±Zj £R, = Vj 


b) 


Vj < real 


±Zj G iR, z'- = Vj 


c) 


Vj, V* complex conjugate pair 


±Zj,±z* G C, 2 pairs = v\*^ 


Table 1. Non-zero complex eigenvalues for real matrix elements 13 = 1. 



Let us now turn to the non-Hermitian generalization. As sketched in 
the introduction we consider a Gaussian two-matrix model and compute 
the spectral properties of the non-Hermitian matrix DC: 



/ 



dCdD exp 



N 



Tt{C^C + D^D) 



(2.5) 



Both C and are rectangular matrices of size Nx{N+v) with real (/3 = 1), 
complex [P = 2) or quaternion real (/? = 4) matrix elements without further 
symmetry. Instead of the N complex eigenvalues Vi of the product matrix 
DC we may again consider the spectrum of the non-Hermitian Dirac matrix 
Ojv C 

feature. The characteristic equation of the real asymmetric matrix DC 



P 



instead. Here the ensemble (/3 = 1) has a special 



remains real. Therefore the solutions of det[A — DC] = Hi^iC-^ ~ 
either real, or occur in complex conjugate pairs. For the non-zero complex 
Dirac eigenvalues Zi we thus get 3 different possibilities listed in Table [TJ 

In the following we will solve a more general non-Hermitian extension of 
the Wishart ensembles than eq. ()2.5p . We add a non-Hermiticity parameter 
/i G [0, 1] that allows to interpolate between the Wishart ensemble eq. ()2.ip 
for /X = and its generalization eq. ()2.5p for /U = 1 which we call maximally 
non-Hermitian. The partition function is defined as 

Zj^^^(/i;{m/}) ~ J dCdD exp -Ni]+Tr{C^C + DD^) (2.6) 



are 



X exp 



-N7]_Ti{DC + CtZ)t)J Yl det[l2N+^mf + p] , 

/=i 

(2.7) 



In the second line we have added a product of characteristic polynomials 
to the weight function, these are called Nf quark flavors in the language 
of QCD. All complex (and real for /3 = 1) eigenvalue correlation functions 
as well as the partition functions are known in the presence of these terms. 
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see [32l [U EJ [3l HI El [9l [M] ■ For simplicity we will mainly restrict ourselves 
to the so-called quenched case Nf = in the following, although Nf ^ 
is particularly interesting for the QCD application as it may lead to a non- 
positive overall weight, the so-called sign problem. For more discussion we 
refer to these references, as well as to [38j. 

The fact that the ensembles eq. (j2.6p indeed extrapolate between the 
maximally non-Hermitian and Hermitian Wishart ensembles can be seen as 
follows. The change of variables 

C = ^+iH^, D = ^^ + iH^'^, (2.8) 



leads to uncoupled Gaussian weights in the matrices $ and ^ (see e.g. 
For = 1 this change of variables is trivial and we are back to eq. (|2.5p . 
Aspects of this simpler model at /x = 1 were also treated in [H [2H [261 112] • 
For /i = 0, D = C\ and so p becomes Hermitian and independent of ^ 
which thus decouples, to give the standard Wishart ensembles eq. ()2.ip in 
terms of matrix <I>. For completeness and later use we also state the elliptic 
extension of the Ginibre ensembles corresponding to eq. ()2.6p (at Nf = 0): 

^NUn^l^)- l^'J exp[-Afa{r?+T:(Jjt)-r?_Tr(j2 + Jt2)|] , (2.9) 

where J ^ is an NxN non-Hermitian matrix. In [20] we have a = {l+v"^) 
and r = {1 - v'^) / {I + v'^) . 

Let us now turn to the complex eigenvalue representation of eq. (|2.6p . 
defining the joint probability distribution 

N 

Here we present the Wishart picture only. For the complex and quaternion 
real case we obtain 

N 



4^;'^(m) = I \{<fz,wf='\z^)\^N{{z])\\ (2.11) 
^V^if^) = lld'z,wif'='Hzj){z,-z*)A2N{{z,z*}), (2.12) 



where we have defined the weight function 



wlP\z)^\z\2'^K^,/2mrj+\z\)exp[Nrj^{z + z*)]. (2.13) 
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The appearance of the Bessel-X function can be understood as follows, just 
considering scalar variables. While the sum of two random variables (the 
real and imaginary part) is again Gaussian - this corresponds to the Ginibre 
case - here we consider the product of two Gaussian random variables, which 
is distributed with respect to Bessel-i^o as one case easily convince oneself. 
For more details of the derivation of the Jacobians we refer to |321 [2], Note 
that for /3 = 4 the Jacobian is different from |AAr({z})|^ as one would 
expect for a standard Dyson-gas For an interpretation in terms of 

charged particles we refer to |18J. Very recently an alternative derivation of 
the Jacobian for /3 = 2 has been presented in |26) . 

In the third and most difficult ensemble (3 = 1 with rectangular real 
matrices one has in principle to sum over all possible combinations of com- 
plex conjugate and real eigenvalues (and purely imaginary ones in the Dirac 
picture) [6j. We only quote the result for the partition function of [8] where 
a factorized form was shown, reading as follows: 

„ 2n „ n 

= j dy"" h^y)l[ J d^ZkllF,{z2,-i,Z2,) A^+2n(y,{4) • 

(2.14) 



The anti-symmetric weight function is defined as 



F,{zuZ2) = i5.(^i,^2)(e(9mzi)-e(9mz2))5'(^2-^i*) (2.15) 
+-hy{zi)hy{z2)5{^m. zi)5{^m. Z2)sgii{^e Z2 - Kezi) . 

It is given in terms of the following two weights for real and complex eigen- 
values respectively: 

hy{x) = 2|xr/2i^^/2(2iVT?+|2;|)exp[27V7/_x] , (2.16) 



°° dt 



g,{zi,Z2) = 2|ziZ2r/'exp[2iVr,_(zi + z2)] / ^e-^^'''i*(^?+^2)-i^ 

Jo t 

{m'^riltziZ2) erfc (2N7]+Vt\z2 - zi 

These two weights are related by 

lim gy{z,z*) = K{^ezf . (2.17) 

While both contain parts of the weight function eq. (j2.13p with /3 = 1 
there, in particular the weight for complex eigenvalues is more complicated 
here. In eq. (|2.14|) valid for both even and odd N = 2n + n ^ 'H the 
notation indexed by x means that for even (odd) N with x = (1) the 
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additional integration over the real eigenvalue y is absent (present). The 
fact that an absolute value around the Vandermonde determinant is absent 
here is due to the ordering enforced by the weight F in terms of the sign- 
and 0-functions. For details on the computation of the Jacobian we refer 
to [61 [Mi- 



3. Correlation functions at finite N: (skew) orthogonal Laguerre 
polynomials in the complex plane 

We now turn to the computation of complex (and real) eigenvalue cor- 
relations functions. The A;-point density correlation functions are defined 
as 

<(.!,..., ^.) - I U . (3.1) 



j=k+i 



The map to the density correlations of Dirac eigenvalues (p) is then given 
by 

k 

Rf^\z,,...,z,) = 2''^ll\z,\'R^S\zl...,zl). (3.2) 

i=i 

For all three ensembles these can be solved in terms of a kernel defined in 
term of Laguerre polynomials in the complex plane as we will show now. 
Analogous results hold for the Ginibre ensembles with /3 = 2, 4 and 1 in 
terms of Hermite polynomials in the complex plane as shown in [20 1 125 1 [T5]. 
respectively. The case /3 = 1 is again special as the /c-point function will 
consist of a sum of all possibilities of real and complex conjugate eigenvalue 
pairs of a total number k (see e.g. [6j). Other correlation functions can be 
defined and computed in these ensembles as well such as gap probabilities, 
and we refer to [5j for more details. 

For /3 = 2 the /c-point correlation functions can be solved in terms of 
the kernel of orthogonal polynomials, in complete analogy to the Hermitian 
case: 



k 



<-')(zi,...,Zfe) =r\wif^='Hzi) det /cS^-')(z„z,) . (3.3) 

1=1 - ■'- 

For the weight eq. (|2.13p it is given through the orthonormalized Laguerre 
polynomials in the complex plane [321 [1] 



fc=0 



(3.4) 
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The Laguerre polynomials satisfy the following general orthogonality rela- 
tion 



/ dz^\z\'^ Kjj {a\z\) eiqi 
Jc 



-iz + z 



Lj(cz) LI {cz*) = hf5j 



for a > 6 > with c = {a? — 6^)/(26). The squared norms read 



7r(j + v)\ (a\'^i 
a jl \b 



2a 



62 



u+l 



(3.5) 



(3.6) 



A short proof for this relation stated in |32j can be found in Proposition 1 
in [7] (see also appendix A of [2j for an earlier proof). The simplest example 
for a correlation function is the spectral density 



wf=^\z)lc'~S-^\z,z*) , (3.7) 

with the kernel from eq. ()3.4p . see Fig. [TJ At maximal non-Hermiticity 
/i = 1 it reduces to an incomplete Bessel-/ function times the Bessel-if 
from the weight [1] 



o(/9=2)/ N 



Ku{N\z\) ^ 



1=0 



tt{1 + u)\11 



N\z\ 



(3.8) 



For more details about correlations functions including Nf characteristic 
polynomials we refer to [Mill]- 

For /3 = 4 the fc-point correlation functions are given by 

k 



R 



(/3=4) 

'N 



Zk) 



\{w\f^='^\zi){z: - Zl) Pfy = l,...,2fc 



1=1 



JC 



(/3=4) 
N 



{Ui,Uj 



where the arguments of the antisymmetric kernel 

N~l 



(3.9) 



/C 



(/3=4) 
N 



(^1,2^2) 



^ /,(4) 
k=0 "fc 



(3.10) 

run through the set of 2k variables {uj} = {zi, zl, . . . , Zk, z^^}, and the 
Pfaffian of an antisymmetric matric of size 2k is given by the square root of 
the determinant of the matrix. The polynomials inside the kernel are given 
by [2] (with (?2fc+i^(^) determined up to a constant times q2k~'^\^)) 



-(2fc + l)! 



N 



2fc+l 



r2i^ 
^2k+l 



Nz 



(3.11) 



N 



2k k 

E 

j=0 



V\^^ k\ (fc + z^)!(2j)! 2. 
r?+; 22^j! (j + i.)! 



Nz 
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with squared norms 



h 



(4) 



87r/(2A; + l)! {2k + 2u + l) 



+ 



2\ik+2v 



They enjoy the following skew orthogonality conditions 

.(4) |J4)> _ /„(4), (4) > _ ,(4) 



(3.12) 



with respect to the following skew-product 

{h\g)p=, ^ [ cPz wf='^\z) (z* - z)[h{z)g{zr - h{zYg{z)] . (3.14) 



,(4)u(4)n 



, 



(3.13) 



Hence they are called skew orthogonal Laguerre polynomials in the complex 
plane. Again we give the simplest example, the spectral density plotted in 

Fig. m 



w 



(/3=4)/ 



\z){z*-z)IC^^-''\z,z*) , (3.15) 

obtained by inserting the kernel eq. (|3.10p . At maximal non-Hermiticity it 
reduces to [2] 



R 



{/3=4), 
N ' 



z\'-K2,{N\z\){z* -z) 



N' 



2v+2 



vr 



(3.16) 



N-l k 

EE 

k=Q j=0 



A:!(A: + z^)!(z2'=+iz*2j _ ^2i^*2fe+i) 
{2k + 2v + l)\{2k + l)!24ij!(j + u)\ ' 



Because of the prefactor the density vanishes identically along the real 
line, or after mapping to Dirac eigenvalues along the real and imaginary axis. 
For results for the correlation functions including Nf pairwise degenerate 
and nondegenerate flavors we refer to [2j and [3], respectively. 

For /3 = 1 we will only give the results for the fc-point correlation func- 
tions for even N = 2n. They are give by [6j 



R 



(/3=1) 
■N 



{zi,...,Zk) =Pf 



fCN{zi,Zj) -GN{Zi,Zj) 
GN{zj,Zi) -WN{Zi,Zj) 



(3.17) 



l<i,j<k 



In addition to the antisymmetric weight eq. (j2.15p we have introduced the 
following two functions of two complex variables 



Gn{zi,Z2) 



(fzlC^^ ^\zi,z)F^{z,Z2) , 



(3.18) 



Wn{zi,Z2) = -F,{zi,Z2)+ I (fz / cPz'F,{zi,z)l6^-^\z,z')F,{z',Z2) 



Allv2 printed on January 20, 2013 



11 



These are given in terms of the integrals of the antisymmetric kernel [3] 



For odd N the corresponding results can be obtained for example by the 
limiting procedure proposed in [16j, or by following [35j. 

The above kernel can be expressed in terms of the following skew or- 
thogonal polynomials [8| 

qg='\z) = +CUzl (3.20) 
ifkth^) = -C!i,^Az)+N~\l + ^,'f{2k){2k + u)Cli,_,{z) + c'C^^,{z), 

giving the skew-orthogonal Laguerre polynomials up to an arbitrary con- 
stant c'. Here we have defined 



f{zi) g{zi) 

f{z2) g{z2) 



N J " \8fi\ 
The corresponding skew-product is defined as 

{f\g)p=i = -{g\f)0=i = j d'^zid'^Z2F^{zi,Z2)det 

(3.22) 

for two functions f{z) and g{z) that are integrable with respect to the 
weight functions contained in Fi,{zi,Z2) in eq. (|2.15p . Our skew-orthogonal 
polynomials defined above then satisfy 

iQ2k\Q^i^)i = il^^+ilQ^ili)! = OVfc,/>0, (3.23) 
where the /i^^^ > are their positive (squared skew) norms 

= 87r(V)(2A;)! {2k + u)l {8^?ri+/Nf^+^+^ . (3.24) 
The kernel is then reading 

n— 1 ^ 

4':2i(^l>-2) = E TTTT U2U'l)<lih-2) - 'Z?il(-2)gil^(^l)), (3.25) 



A:=0 
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leading to eq. ()3.19p . 

As already mentioned, here the /c-point function in eq. (jS.ip is the sum 
of all possible combinations of complex and real eigenvalues, starting from 
only complex to only real eigenvalues. For example for the spectral density 
plotted in Fig. [3] we have [6] 

R^j^\z)= [ <fz'lC^=\z\z)F{z,z')^R^i,\{z) + 5{y)R^},\{x) . (3.26) 
Jc 

Here we explicitly give the result obtained after inserting the kernel (j3.19p 
and integrating with respect to the two different weights eq. (j2.16p . Denot- 
ing z = X + iy we get in terms of the kernel from eq. (|3.19p 

i?«^(z) = -2i{N\z\Ye^''^-^sgn{y)IC^j^\z,z*)2 r '^e-^NSlt{x^-y')-h 

Jo t 

xK^ (2N'^y]\t{x^ + y2)^erfc (2Ni^+Vt\y\^ , (3.27) 

/oo 
dx'sgn{x - x') |xx'|§e^''-(^+^')K|(iVr/+|x|) 
-oo 

xK^_{Ni^+\x'\)ld^\x,x'). (3.28) 
At maximal non-Hermiticity the kernel eq. (|3.19p considerably simplifies [3] 

Finally in the map to the Dirac picture one has to distinguish between 
substituting complex or real eigenvalues. For the densities we have 

<^)(z) = A\z\^d;^]^{z% (3.30) 
^N,R (■^) ~ 2|x|i?jY^'']g(x^) , (3.31) 
and analogous relations hold for higher /c-point functions. 

4. The large- AT limit: 3 classes of complex Bessel kernels 

In this section we will review aspects of the large-A^ limit for non- 
Hermitian Wishart random matrices. Because the eigenvalues live in the 
complex plane we have more possibilities in taking the large-A^ limit, de- 
pending on the location in the spectrum. In order to illustrate where such 
regions are located in the complex plane we first display the so-called macro- 
scopic (or global) limit of the spectral density. As one can see from plotting 
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the density at finite- for all three j3 = 1,2,4 in the Dirac picture in Figs. 
[T]to[3l the global spectral density follows the elliptic law [22j as their Gini- 
bre counterparts. On top of being constant on an ellipse, various inner and 
outer edges exist when zooming into particular regions. This is in contrast 
to the real density of Wishart matrices where one can only distinguish the 
(soft) edge, the bulk and the origin (or hard edge) region. 

In addition to the location of the spectrum two different large- lim- 
its have to be distinguished for complex eigenvalues, that of strong and 
weak non-Hermiticity [H]. While in the former the eigenvalues fill a two- 
dimensional support, in the latter they only locally extend into the complex 
plane. The weakly non-Hermitian limit turns out to be a one-parameter 
deformation both of the Hermitian and the strongly non-Hermitian limit. 
This makes it an ideal object to study the question of universality. 

In the following we will focus on the microscopic limit at the origin, 
as this is a special feature of Wishart ensembles that their corresponding 
Ginibre counterparts cannot offer. 

We begin once more by recalling the results for the macroscopic spec- 
tral density in Hermitian Wishart ensembles. Here two cases have to be 
distinguished. Taking the large- limit with u = 0{N), such that c = 
N/{N + ly) < 1 one obtains the Marchenko-Pastur distribution [31] for the 
positive Wishart eigenvalues for all three values of /3 



p(P,MP)^X) = -1- V(A - cX.){cX+ - A) , X± = Vc-V2±i (4.1) 



on A E [cX^,cX^] and zero outside. For asymptotically square matrices 
with 1/ = 0(1) at large- (leading to c = 1) it reduces to 



which is just the semi-circle after mapping to the Dirac picture. While in 
the Wishart picture we have a square root singularity at the origin, the 
density in the Dirac picture is regular at the origin. 

The same phenomenon happens in the complex plane. The semi-circle 
density becomes replaced by the elliptic law, a constant density of an ellipse 
[22j . which we state for the elliptic Ginibre ensemble eq. (j2.9p (see e.g. [20j ) 
in variables z = x + iy 




(4.2) 





for Ul + v'f + f^{l + v'f<l 
4 Av^ 



(4.3) 



P 



and zero outside. Away from the origin the mean spectral density of our 
three non-Hermitian Wishart ensembles is constant on an ellipse in the Dirac 
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picture, just as its non-chiral counter parts. This can be seen in the figures 
below Fig. [U - O where we plot the corresponding spectral densities after 
mapping to the Dirac picture according to eq. (|3.30p (and eq. ()3.3ip for 
/? = !)• 




Fig. 1. The Dirac spectral density ^' ^\z) for /3 = 2 from eq. (P?71) for = 10, 
PL = 0.7 and v = Q (left), and v — \ (right). The extra repulsion from the origin 
through the zero-eigenvalue at = 1 is clearly visible. In the large-A limit the 
detailed structure at the origin will only be visible on a microscopic scale. 




Fig. 2. The Dirac spectral density ^ (z) for /3 = 4 from eq. (j3.15p for 
N = 20, /i = 0.7 and v ^ Q (left) and v ^ \ (right). In order to see a smoother 
plateau we have increased N here. The extra repulsion from the origin for v — 1 
is less pronounced here compared to P — 2. 



When mapping back to the Wishart picture using eq. ()3.30p we obtain 
a linear decay inside the support, as was observed in [26} [T2]: 

p^P^^\z)r. const. pW(z)^^^. (4.4) 
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Fig. 3. For (3 — 1 we show two Dirac spectral densities: Rj^ ^ [z) for the 

complex eigenvalues, and (ly^\z — x) {{iy)) for the real (or purely imaginary) 
eigenvalues from eqs. p.27p and (j3.28l) . The parameter values are the same as 
for /3 = 4 in Fig. [2] It is clearly visible that the density of complex eigenvalues is 
repelled from the real and imaginary axis as for /3 = 4, with an apparently different 
profile though. 



Likewise the behavior at the origin, that was found to be constant for > 
and logarithmically divergent ~ log \ z\ for = in the Wishart picture |26| . 
gets bent to zero in the Dirac picture using the map eq. (|c{.3U|) . This can 
bee seen in the figures for all three (3 and the values of = 0, 1. 

In addition, in [Ml 112] the limit v = 0{N), such that q = v/N > 
finite, was performed for /3 = 2 at /i = 1, leading to the density 

on the support, which we give in the Wishart picture here. This is the gen- 
eralization of the Marchenko-Pastur distribution eq. ()4.ip into the complex 
plane. 

Let us now turn to the microscopic large- limit. As can be seen from 
the figures several regions could be investigated by zooming into the outer 
edges or inner edges along the axes. However, in the following we will 
restrict ourselves to the microscopic origin limit, as this limit is particular 
to the Wishart ensembles. 

We begin with the microscopic origin limit at strong non-Hermiticity 
{S). Because strong and maximal non-Hermiticity at = 1 are related by a 
simple rescaling of the complex eigenvalues, z — )■ 2rjj^z (see e.g. [201 E]), we 
restrict ourselves only to the latter. This limit is defined such that the non- 
Hermiticity parameter ^ is not scaled with N , only the complex eigenvalues 
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get rescaled: 

pf'f'Xi) = lim li?!^^'^)(z = e/\/iv) , lim ^z = i fixed, (4.6) 

and all higher order correlation functions are rescaled accordingly. 
For /? = 2 we obtain easily from eq. ()3.8p [Ij 

pf^\0 = -\C\'KAmm\')^'^-. (4.7) 
vr vr 

It only depends on the rescaled modulus |^| and is thus rotationally in- 
variant. We also give the asymptotic value of the microscopic density that 
matches the constant value of the macroscopic density in the Dirac picture. 
For /3 = 4 we obtain from eq. (|3.16p after a non-trivial calculation [2] 

pf^ho = ^-^ic'-e)K2u{m (4.8) 



1 hA^iej ^^^u I vi 

2 



Compared to /3 = 2 the density is no longer rotationally invariant, reflecting 
the repulsion of the complex eigenvalues from real and imaginary axes. 

For /3 = 1 we obtain two densities from eqs. ()3.27p and ()3.28p . for the 
complex eigenvalues and the eigenvalues on the real and imaginary axis [6]: 

Pcfii) = ^m{'^He)){-^){e -e')l\eiu[m'') 

^ j^j e-«'+«")*-i* {2\i\H) erfc {2Vt\^m{e)\) , (4.9) 
Pml^O = \\i\K^^i\i\'')[j^ dx'\e-x'\ K.^{\x'\) h{2i^) 

+ f dx' \e - x'\ K^{\x'\) M2Cy^\)) , (4.10) 

where in the second equation ^ G M (or ^ G which is identical here). 
Once again the density in the complex plane is not rotationally invariant, 
and vanishes along the real and imaginary axis. 

We now turn to the microscopic origin limit as weak non-Hermiticity 
(W) first introduced for the Ginibre ensembles eq. (|2.9p in [19j. It is 
defined by both rescaling fi with A^, as well as the complex eigenvalues 
(with a different power compared to the strong limit above): 

pf'^\0 ^ Jim ^<^)(. = ^/iV) , (4.11) 
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with 



lim 2iVV^ 

N^oo 



a 



(4.12) 



fixed. All higher order correlation functions are rescaled accordingly. 
For /3 = 2 we get [32] 



" 27ra2l^l'^'^(4a2 

5(9me)^(J.(e)'- J.-i(OJ.+i(0) • 



(4.13) 



It agrees with the density obtained from a non-Hermitian one-matrix model 
[26j. In the second line we indicate that in the Hermitian limit a — >• the 
density reduces to the known universal Bessel density times a delta function 
in the imaginary part. In the opposite limit a — t- oo the density at strong 
non-Hermiticity given in eq. (j4.7|) can be recovered, and we refer to |T] for 
more details. 

For /3 = 4 the weak non-Hermiticity limit yields [2] 



Pw 



32a4 



2u 



2a2 ) 



ds j — e 
Jo Vt 



■ (4.14) 



Once again in the Hermitian limit a — )• the prefactor in front of the integral 
reduces to a delta function in 9m(i^), and the integral matches the known 
universal Bessel kernel for /? = 4 (which also has a simpler representation 
expressed in terms of the (3 = 2 density plus a single integral) . 

Finally we give the weak origin limit for /? = 1 which was obtained most 
recently. The complication arises from the fact that for the real eigenvalues 
both at // = in the real Wishart ensemble, and at /i 7^ for our non- 
Hermitian extension, the large- limit of the kernel and the integration in 
eq. ()3.28p do not commute. For more details we refer to [9l [3l] and only 
quote the answer valid for ^ G M and ^ G 



[sgn^2] 2 
1 

'320F 



2 



2dy 



V2a 



MO + 



[sgn^2 
2{V2a) 



V /.I 



J,+i(se) - 2c? s (E+is) - E.is)) MsO 



(4.15) 
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where we have defined the kernel in the weak Hmit 



Jo 



297ra2 




(4.16) 



the rescaled real weight function 



V — ■ 

2 Vl6a2 



(4.17) 



and the exponential integral 




oo 



(4.18) 



The weak density of complex eigenvalues presents no such difficulties [6]: 



In the Hermitian limit a — >• the density of complex eigenvalues as well 
as the density along the imaginary axis vanishes. Only the density of real 
eigenvalues will build up the known density of real eigenvalues in the Her- 
mitian Wishart ensemble for /3 = 1 [6]. 



In this short article we have reviewed the solution of three non-Hermitian 
extensions of Wishart ensembles of random matrices with real, complex 
or quaternion real elements. In all three cases the eigenvalue correlation 
functions are expressed in terms of the kernel of (skew) orthogonal La- 
guerre polynomials in the complex plane, depending on a non-Hermiticity 
parameter that allows to interpolate between Wishart and maximally non- 
Hermitian Wishart ensembles, for finite and infinite A^. At the origin these 
ensembles are very well understood, and we gave the corresponding three 
Bessel kernels in the complex plane at strong and weak non-Hermiticity. 
At the various inner and outer edges we expect that the behavior of the 
corresponding Ginibre ensembles (see e.g. in [27]), and the ensembles we 
considered here is be the same and thus universal. This was shown for ex- 
ample in [7| for the generalized Airy kernel at the soft edge along the real 
line, or in in the rotationally invariant case (regime II (iii)) at the outer 
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edge, both for (3 = 2. The latter behavior was also found for the outer edge 
density of random contractions [27]. 

At present a deeper guiding principle to prove universality for non- 
Hermitian random matrices with a larger class of weight functions, such 
as quasiharmonic potentials [39j is lacking. A first step towards universality 
could be to show that different Gaussian models lead to the same answer as 
mentioned above, or as it was shown in yet another example for the Bessel 
kernel in the weak limit at /3 = 2 eq. (|4.13|) . starting from a Gaussian one- 
or two-matrix model in [36J and [32J, respectively. We hope that these first 
steps will ultimately lead to a deeper understanding of the issue of univer- 
sality. 

Acknowledgments I would like to thank all my coworkers for collab- 
orations on this subject, as well as the organizers of this workshop for the 
very stimulating atmosphere. 
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